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We have searched 4.5 square degrees of archival HST/ACS images for cosmic strings, 
identifying close pairs of similar, faint galaxies and selecting groups whose alignment is 
consistent with gravitational lensing by a long, straight string. We find no evidence for 
cosmic strings in five large-area HST treasury surveys (covering a total of 2.22 square de- 
grees), or in any of 346 multi-filter guest observer images (1.18 square degrees). Assuming 
that simulations accurately predict the number of cosmic strings in the universe, this non- 
detection allows us to place upper limits on the unitless Universal cosmic string tension of 
Gfi/c^ < 2.3 X 10~^, and cosmic string density of Vts < 2.1 x 10~^ at the 95% confidence 
level (marginalising over the other parameter in each case). We find four dubious cosmic 
string candidates in 318 single filter guest observer images (1.08 square degrees), which we 
are unable to conclusively eliminate with existing data. The confirmation of any one of these 
candidates as cosmic strings would imply Gji/c^ ^ 10~^ and ~ 10~^. However, we es- 
timate that there is at least a 92% chance that these string candidates are random alignments 
of galaxies. If we assume that these candidates are indeed false detections, our final limits 
on Gji/c^ and Q^s fall to 6.5 x 10"'' and 7.3 x 10~^. Due to the extensive sky coverage of 
the HST/ACS image archive, the above limits are universal. They are quite sensitive to the 
number of fields being searched, and could be further reduced by more than a factor of two 
using forthcoming HST data. 
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1 INTRODUCTION 

Cosmic strings, which are topological defects formed in the very 
early universe and persisting to the present epoch, are a quite 
generic prediction of many modern cosmological models (Kibble 
1976; Vilenkin & Shellard 1994; Hindmarsh & Kibble 1995). As 
microscopic objects, their internal structure can be characterised 
by a single number, the string tension fic^ which we express in 
unitless form Gji/c? . Numerical simulations of this phenomenon 
have converged on a picture where 10 horizon- scale strings are 
consistently predicted to lie within the observable universe (see 
Polchinski 2007, for a review). This suggests the possibility of a 
direct detection. 



Currently, the best limits on the properties of a putative net- 
work of long strings come indirectly from studies of the CMB 
power spectrum and pulsar timings. The CMB power spectrum is 
now sufficiently well-modeled so as to permit only a very small 
fraction of the power in temperature fluctuations at the last scat- 
tering surface to be due to strings of any type (Pogosian et al. 
2004), and Pogosian et al. (2006) sets 95% condfidence limits on 
of Gii/c^ < 2.7 X 10~^. Pulsar timing experiments like those 
described in Kaspi et al. (1994); Damour & Vilenkin (2005) give 
tighter limits, but make model-dependent assumptions on the rel- 
ative densities of string loops and the long strings we search for 
here. Additionally, pulsar timing limits rely on string oscillations 
to produce gravitational radiation in observable frequency ranges 
(Accetta & Krauss 1989; Mack et al. 2007). But some cosmic 



* E-mail:ericeric @ slac.stanford.edu 



2 Morganson et al. 



string models decay via Goldstone bosons or other field couplings 
(Polchinski 2004b). 

Direct searches for cosmic strings provide an important com- 
plement to those statistical studies. Direct searches in the CMB rely 
on the movement of strings, which introduce a very small appar- 
ent differential between light from other side of the string (Kibble 
2004): 

— — = SttG/x/c smz — - (1) 

z c Ds 

Where sin i is the projection of the string along the line of 
sight, Ds is the angular distance from the observer to the CMB, 
and Dds is the angular distance from the string to the CMB. This 
redshift is potentially observable using the CMB as a backlight. 
The instrumentation required to make these observations at high 
resolution is still in its infancy (Zwart et al. 2008; Muchovej et al. 
2007; Stewart & Brandenberger 2009; Fowler et al. 2007; Bouchet 
2007). Previous attempts at lower resolution have not resulted in 
any detections, and produced upper limits on the string tension of 
3.7x10"^ (Sazhina et al. 2008; Lo & Wright 2005; Jeong & Smoot 
2007). This induced redshift is of order 10~^°, too small to be de- 
tected in the optical spectroscopic galaxy surveys, where it would 
be dominated by galaxy motions. 

Cosmic strings are also directly observable in optical imaging 
surveys via their gravitational lensing effect (e.g. Vilenkin 1984; 
Sazhin & Khlopov 1989). Cosmic strings produce a unique lensing 
signal: two identical images of a source separated by 0: 

= S7vGfi/c^smi^, (2) 

Where Ds is the distance between the observer and the source and 
Dds is the distance between the string and the source. 

In our first paper (Gasparini et al. 2008, Paper I), we showed 
how high resolution optical imaging surveys are capable of produc- 
ing direct limits on the cosmic string tension competitive with the 
indirect detection from the CMB power spectrum analysis. For the 
tensions allowed by the CMB power spectrum analysis, this sepa- 
ration is < 1", and the cross-section for lensing is similarly small. 
This realisation immediately drives us towards the need to use the 
most numerous distant objects — faint galaxies — as backlights, and 
to focus on surveys with high angular resolution. 

Early optical cosmic string searches focused on large image 
separations and bright galaxies, but no candidate string lenses have 
survived inspection at higher resolution (Agol et al. 2006; Sazhin 
et al. 2007). Ground-based optical surveys are hampered by limited 
resolution and cannot probe Gfi/c? < 10~^ (de Laix 1997). Con- 
versely, the recent searches by Christiansen et al. (2008) using the 
Hubble Space Telescope (HST) GOODS survey images and Timo- 
thy et al. (2009) using HST COSMOS survey images did have the 
angular resolution to show that no string with G/n/c^ > 3.0 x 10~^ 
crossed their survey fields, but did not have the sky coverage needed 
to be able to place any universal limit on the string tension. 

We need both large survey area and high resolution to set uni- 
versal limits on string tension. For the case of HST, our best cur- 
rent source of high resolution imaging survey data, the solid angle 
needed to give a universal limit on the string tension is a few square 
degrees, provided this area is cut up and distributed evenly across 
the sky. As discussed in Paper I, a contiguous survey of compara- 
ble area would provide much weaker limits due to the clustering of 
cosmic string events around the roughly 10 long strings in the sky. 

In this work we search the HST/ACS archive, whose images 
match all three of the angular resolution, area and pointing dis- 



tribution criteria, and so produce the strongest direct gravitational 
lensing limit on the cosmic string tension and density. 

This paper is organised as follows: in Section 2, we briefly re- 
view the relevant basic theory of cosmic strings. We then discuss 
how searching for string lensing events can be used to set new lim- 
its (Section 3). Regarding our actual search for string lenses, we 
explain our use of simulated string lensing images in Section 4, 
and in Section 5 we describe our string-finding technique. After a 
brief description of the archival dataset in Section 6, we present our 
results in Section 7 and the corresponding limits on the string ten- 
sion and density in Section 8. Finally, we discuss the implications 
of this work in Section 10. 



2 STRING LENSING THEORY 

When searching for cosmic strings via their gravitational lensing 
effect, we only need to understand a few basic facts about them. 
We must know what type of multiple image system a cosmic string 
produces, and we must understand what string evolution tells us 
about the number of strings in the local (zd ^ 2) universe. 

Cosmic strings, whether topological defects or stretched fun- 
damental strings, can produce two identical images separated by 
in Eq. 2 so long as they do not "cut" the source they are lensing and 
so long as they are straight on scales comparable to the image sepa- 
ration. By "cutting" a source, we mean that the redshift-dependent 
strip of sky which the string copies only includes a part of a source, 
and the source in incompletely copied. We study this effect via sim- 
ulation in Section 4.2. We discuss our assumption of straight strings 
below. 



2.1 String network topology: straight strings and loops 

We can only observe a string by its gravitational lensing effect if 
it crosses one of our search fields; we then expect it to be more 
readily detectable if it is straight on the scale of the field. We 
consult the theory and simulations of cosmic strings to determine 
whether these are reasonable possibilities. Horizon scale strings 
are a generic prediction of string network simulations (Polchinski 
2004a), but we must consider the total angular length of observable 
string and, and to what extent small scale structure (kinks and cur- 
vature) will affect the image pair alignment. We also explore the 
possibility of detecting string loops. 

Simulations (Allen & Shellard 1990) tend to agree that during 
the matter dominated epoch and beyond, the long string density is: 

p, = ylfGp/c^pm; ^ = 60 =b 15, (3) 

where ^ represents the average length of long, straight string 
within a cosmological volume, the concentration of string. These 
simulations assume that strings "intercommute" with probability 
of 1. Lowering this probability raises the simulated value of ^, but 
we assume unity intercommutation probability, because it is theo- 
retically favoured. 

We can relate ^ to the number of cosmic strings within some 

redshift limit. Substituting in pm = ^^r the co-moving 

matter density, we have: 

^ = *n^M (4) 
We divide by 2Dd (where Dd is the angular diameter distance to 
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the string) and obtain the concentration of string in angular units, 
accounting for projection effects: 



< 6>. > 
dV 



(5) 



Multiplying by the cosmological volume element and integrating 
over redshift yields the total expected angular length of string on 
the sky, in radians, < Os >: 



<0s> 

<0s> 
<0s> 



dz 



^0 

Jo 4c^nm{l + Z)^ + 

0.35^, for z^ax = 2 



167rc2D(z) dz 



dz 



(6) 
(7) 



Os ^ TV corresponds to a single straight string stretching across 
the visible sky. Setting ^ = 60 gives us 21 radians of string within 
redshift 2. This in turn corresponds to roughly 7 straight strings 
with redshift < 2. Even survey fields as large as a few square de- 
grees would most likely not be crossed by a string, and would there- 
fore have no chance of seeing a long string regardless of survey 
depth or resolution (see Section 3.1 for details). A survey of the 
same area comprising a large number of smaller, well-distributed 
fields stands a much better chance of "hitting" a string. 

Even if a string crosses one of our search fields, detecting it 
may be nontrivial. Cosmic strings that are straight on scales com- 
parable to the field size should be readily identifiable: groups of 
aligned galaxy pairs make good candidates for being multiple im- 
age systems generated by the string. Fields of ~ 1 arcmin would 
then result in the requirement that, for it to be detectable, a cosmic 
string be straight on Mpc scales. 

What can we infer from the string network simulations about 
the likely straightness of horizon-scale strings? Suppose that cos- 
mic strings follow random walks. The total length of a string within 
the horizon of a string undergoing a random walk is roughly d^/L 
where dh is the horizon scale and L is the characteristic length 
of the random walk. Assuming a string tension Gji/c? > 10~^, 
the random walk length, and similarly the radius of curvature of 
the string, must be of order the dh or ^ would quickly exceed 60 
and ps would be large enough to affect the CMB spectrum (see 
Section 3). So L must be large and this random walk behaviour 
only affects the Mpc scale perturbatively. The simulations do show 
that a string's self-interactions and its interactions with other strings 
can induce smaller scale "cusps" which propagate along the string. 
However, at least some models suggest that these cusps will leave 
strings that are essentially straight on Mpc scales (Rocha 2008). 

String evolution networks also produce loops. Loops smaller 
than a few Mpc would not be detectable as simple alignments of 
image pairs. Mack et al. (2007) provides a framework for detect- 
ing string loops with proposed high resolution radio surveys. Many 
models and simulations set the scale of cosmic loops (Caldwell & 
Allen 1992) to be of order the horizon scale, and so it is possible 
that many cosmic string loops are straight on the Mpc scale, but we 
do not include them when marginalising over string concentration 
in our calculations. 



2.2 String lens observability 

In this section we summarise the various observable features of cos- 
mic string lenses, and thus motivate our particular search strategy. 

We might hope to detect a cosmic string using a single, obvi- 
ous string lensing event, where the images are demonstrably trans- 



lated copies of each other (see Agol et al. 2006, for such an anal- 
ysis). It is only possible to determine if a single pair is the result 
of string lensing if the pair separation is several times larger than 
the PSF, and if the sources are bright and well-resolved. The large 
separation requirement would directly limit our ability to probe the 
small (subarcsecond) string tensions of interest. The bright source 
requirement, and the fact that there are exponentially fewer bright 
sources than dim sources, would limit our ability to probe small ^. 

The key observable feature of a cosmic string is therefore that 
it produces many pairs of images, as it cuts across many back- 
ground sources (Huterer & Vachaspati 2003; Oguri & Takahashi 
2005). We could look for local overdensities of well-matched (in 
magnitude, ellipticity and orientation) pairs, as is expected to occur 
when a kinked or coiled string lies in the field. This is the approach 
taken by Christiansen et al. (2008), and only works when the num- 
ber of string pairs (cx: ^G/i/c^) is large enough that string lensing 
can be detected against the very high background of pairs described 
by the small-scale correlation function. 

However, on arcsecond and sub-arcsecond scales the back- 
ground sources are highly correlated (Morganson & Blandford 
2008). The increase in the abundance of apparent galaxy pairs due 
to the presence of a cosmic string is, as we shall see, quite small. 
Only for the expected straight strings from the previous section do 
we expect a measurable phenomenon: the appearance of groups of 
apparent galaxy pairs aligned with each other. 

The similarity of the members of each aligned pair should still 
be a useful indication of string lensing. At image separations of 
0.3^^ we are probing down to the scale of the source size itself, 
and so many multiple image systems will contain incompletely 
copied sources. If these images were well-resolved, we could de- 
tect a sharp edge, but small, incompletely copied sources look un- 
remarkable. Pairs with incomplete images can, in rare cases, have 
significant differences in their magnitudes and ellipticities. In Sec- 
tion 4.2, we quantify these differences. 



3 ESTIMATING THE EXPECTED NUMBER OF STRING 
DETECTION 

In order to set limits on string parameters, we calculate the expected 
number of string detections in a given dataset, A^string((^/i/c^, 

Each (assumed rectangular) field is defined by three parame- 
ters, two spatial extents Oi and O2, and its limiting magnitude, m. A 
given string is defined by its tension, Giijc? , its inclination along 
the line of sight, i, and its redshift, z^. The intersection between the 
string and the field is defined by ^c, the projected angular length of 
string that crosses the field. 

To calculate A^string, we must first calculate the probability 
that a randomly oriented string crosses a field with crossing length 
Oc-> -Pcross (^1, ^2, ^, sin z, 2;d, ^c). We then calculate the probabil- 
ity that a string which crosses our field with length Qc is de- 
tected, Pdetect (vfi^ Gfi/c^ , sin i, Zd, Oc). Knowing these two quan- 
tities, we can integrate over all possible strings to obtain A^string for 
a randomly oriented field: 



string 



/r /-TT/j 
dOc / dzd / sin^ 



■ cross -^detect 



(8) 



= .Pc 



cross -^detect 



In Section 3.3 below, we discuss how we transform our expression 
for the probability of a detection in a single field to our 95% detec- 
tion limits for a multi-field survey. 



4 Morganson et al. 



3.1 The Probability of a String Crossing a Field 



We define Pa 



ftross(^i,^2,^,sini,2;d,^c) as the proba- 



bility that a string at redshift Zd will cross our Oi x O2 field with 
length 6c at inclination angle, z, to the line of sight. We assume that 
our fields are much smaller than a radian and that our strings are 
straight on scales longer than a few arcminutes (the size of a sin- 
gle field). We break down the probability of a string crossing our 
field with overlap length Oc into two terms: the probability of the 
field crossing a randomly oriented line that is 1 radian long and the 
projected length of string at redshift Zd (in radians): 



-Prad {Ol ,02,0c 



Ps sinz D^dD 



dz 



(9) 



Prad is a quantity derived solely from two dimensional spherical 
geometry without any cosmological input. It is the probability of a 
randomly positioned Oi x O2 rectangle being crossed by a randomly 
oriented 1 rad arc with an overlap length of Oc. The angular length 
1 radian is arbitrary and chosen for mathematical convenience. We 
provide its (somewhat cumbersome) geometrical derivation in the 
appendix. 

The remaining terms represent the expected projected angular 
length of string on the sky at every redshift and the cosmological 
volume element. If we plug in our model value of ps (equation 3) 
and the Acdm volume element, we obtain: 

Pcross - PMO,M)^ 2C V^^(l+.)3+n, ^''^ 



3.2 The Probability of a String Crossing being observable 

Assuming that a string crosses a field with overlap 0^ we define 
the probability of that string producing a detectable lensing signa- 
ture as Pdetect = Pdetect {m, Gp/c^ siiii.Oc, Zd). To obtain this 
function, we must calculate the expected number events for our 
string and field. This is just a (zs -dependent) lensed area times a 
(2;s-dependent) source density, integrated over Zs. 

The lensing cross sectional area cr of a string is its length, 0^ 
times the string lensing width in equation 2: 

(7 = (T{zd, Zs, Gp/c^ ,smi, Oc) = OcSnGp/c^ sinz (11) 

We model the source density, n{zs,m), using the redshift distri- 
bution distribution from Leauthaud et al. (2007) and a fit to the 
magnitude distribution of COSMOS: 



n{zs,m) = 

no{m) = 
zo — 



no{m)—^ey^oJ 

OZq 

0.67(m-34.9) -2 

e ^ ^ arcsec 

0.13m - 2.3 



(12) 



This equation assumes magnitudes m observed in the F814W filter. 
For fields imaged in a different filter, we measure the source den- 
sity and assume the F814W limiting magnitude that would yield an 
identical density. 

We multiply the cross sectional area by the (redshift- 
dependent) density of sources and integrate over source redshift 
to obtain the expected number of lensing events 



I dzs 

Jz^in Jo 



dm dzscrn{zs,m) (13) 



Here, ^^min is the minimum source redshift that will produce a re- 



solvable lensing event: 

Ds (^Zmin) 



Dds {^d, Zmin) 



= S7V 



Gp sin i 

i-'res 



(14) 



^res is the "effective resolution," or the angular separation at which 
actual faint pairs in our survey can realistically be deblended. Mor- 
ganson & Blandford (2008) found ^res ~ 0.3", and we use this 
value when deriving detection limits. 

Most of our events will involve faint sources that are fairly 
weakly clustered on arcminute scales, so we can assume a Poisson 
event rate. We require three aligned pairs of similar galaxy images 
to claim a string detection, so we calculate the probability of mak- 
ing three detections or more: 

p, _ 1 „-<A^events> A i ^ Ar ^ , < iVevents >^ ^ 
-^detect — 1 — e llH-< iVevents > H ^ I 

(15) 

When probing a potential string detection, we incorporate the 
information we get from the detection into the above analysis. The 
expected number of lensing events becomes 

/-^max r^lim 
< Nevents >= Oc / dZs / dm dZsCTTl^Zs, TTl) (16) 

Where, Zmin and Zmax are the observed source redshifts that will 
produce the minimum and maximum observed lensing separations: 



Ds (^Zuiin) 
Dds(^Zd, Zuiin) 
^max J 



= Stt 



Gp sin i 
Gp sin i 



(17) 



^max J ^ 

In addition, we define Pobservation to include the observed number 
of detected pairs: 

-<A^events> < -^events >-^--^°°— ' 



Pobservation — 6 



iVol 



(18) 



3.3 Expected number of string detections and string limits in 
a multifield survey 

In this paper we aim to set limits on Gp/(? using multi- 
field surveys, and as an intermediate step, we must calculate 
A^string survey (C/i/c^, ^), thc cxpcctcd uumbcr of String detec- 
tions across our complete dataset. Nstring{Gp/c^ ,'^) must be 
calculated numerically for every field at each tension Gp/c^, 
but it is linearly proportional to ^. For each field, we calculate 

-^string 



(Gp/c^, ^ = 1), and then compute Nstring{Gp/c^ , ^) 



strins{Gp/c',^ = 1). 

The expected number of string detections of multiple fields 
will add linearly so long as the fields are randomly distributed, so 
we have: 

-^string survey (G^^/c^*) = Yl ^string i(GM/c^*) (19) 

fields i 

= ^ A^string survey (Gp/c , ^ = 1) 

To see what Gp/(? and ^ give us 95% chance of detecting a string, 
we assume a Poisson distribution for the number of detections. A 
random variable taken from a Poisson distribution with a mean of 
log 20 = 3.00 only has a 5% chance of being 0. So for a given 
Gp/c^ , we set: 



A^. 



95% 



string survey 



95% N 



log 20 



95% 



log 20 



A^string survey — 1) 



(20) 
(21) 
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To compute proper Bayesian confidence limits, we must must 
assume prior distributions for Giijc? and ^. We assume Gfi/c^ is 
log-uniformly distributed between 10~^ and 10~^, and we take a 
logarithmic version the ^ distribution from Eq. 3: 



P{L) = -^e ^ (22) 
V27rcr 

L = log(^) 

Lo = log (60) 

. _ log(2) 



The probability of a null detection is e~^^*'-^^s^^^/''^'*^ 
and the probability of observing a detection in 
one field and no detections in all other fields is 

e-^^*^i-s(^^/"''*)Pobservation(GM/c^^). Wc USC the aboVC 

priors, these formulae and a standard Bayesian formulation to 
produce confidence limits in Section 8. 



4 SIMULATED IMAGES 

As we discussed in Paper I and Section 2 above, to detect a cos- 
mic string lensing event one must be able to detect close pairs of 
faint objects. Developing and testing search algorithms, and then 
characterising their respective selection functions, requires sets of 
simulated images in which we know the location of the sources 
and the strings. These mock images are based on those used in 
the faint source correlation function measurement by Morganson & 
Blandford (2008). These simulations reproduce the magnitude, area 
and ellipticity distributions of the faintest galaxies in HST/ACS 
GOODS images; we adapt them for cosmological settings with cos- 
mic strings by including magnitude-dependent clustering, stochas- 
tic redshift assignment and a string lensing effect. We provide a 
very brief outline of the image simulation procedure we use in this 
section, and refer the reader to the paper by Morganson & Bland- 
ford (2008) for more details. 



4.1 Image Production 

Our standard simulated images are square, consisting of 8192 x 
8192 pixels. Each pixel is a 0.03" square, so our mock frames are 
246" on a side. In each mock image we place a single straight string 
(as described in Section 4.1), with known orientation. We use these 
images to test our string detection algorithms and to measure the 
detection probability. We also produce a set of narrow 10" x 246" 
images, each containing a single string, in order to produce a high 
density of multiply-imaged faint galaxies and allow us to quantify 
the string lensing effect on the galaxy images themselves (see Sec- 
tion 4.2). 

We choose source redshifts and magnitudes using the distribu- 
tions in equation 12. Each source is assigned an independent red- 
shift, based solely on its magnitude, and printed into the image with 
the nearest discrete Zs = (0.1, 0.3, 0.6, 1, 1.5, 2, 2.5, 3, 3.5, 4, 5, 6). 
We have one string per image. The strings are located at a constant 
redshift (the analog of the thin lens approximation). Each Zs image 
is lensed by copying and shifting half of the image by the amount 
given equation 2. All simulated strings are absolutely straight and 
vertical. We add the images at each Zs to produce an image with all 
sources. 

Finally, we convolve the lensed and combined images with 



a point spread function (PSF) determined by an average of many 
model PSFs calculated using the "Tiny Tim" software package. 

We then add a noise background which has been convolved 
with a separate PSF that mimics the image combination process. 
We vary the noise amplitude to produce images of various desired 
limiting magnitudes. 

In Fig. 1, we show 5 examples of typical string-lensed image 
pairs from one of our simulated fields. Overlaid are ellipses repre- 
senting the orientation and ellipticity measured for each source us- 
ing the SEXTRACTOR program (Bertin & Arnouts 1996). We see 
that some pairs will have well-measured and well-matched mag- 
nitude, size, ellipticity and orientation, and that these should be 
closely-matched between conjugate images in the multiple-image 
system. However, due to the prevalence of faint, poorly-measured 
sources, and to strings splitting sources, most pairs are not as read- 
ily recognizable as string lensed pairs. We quantify this below. 



4.2 Properties of simulated string-lensed objects 

One important application of our simulations is to compare the sta- 
tistical properties of string lens multiple-images with those of un- 
lensed, background galaxy pairs. An idealised string lensing event 
would produce two identical galaxy images, but incomplete copy- 
ing of sources leads to non-ideal string-lensing. Still, as we show in 
this section, string-lensed pairs tend to be more similar than random 
pairs. 

We compare catalogs of string-lensed image pairs and un- 
lensed close galaxy pairs generated from a large number of sim- 
ulated images. These images are the long, narrow strips mentioned 
in Section 4. We made a set of thin (10" x 246") simulated images 
with a string running along the long axis. For each combination of 
lens redshift za = (0.1, 0.5, 1, 1.5, 2, 3), SnG/n/c^ =(0.3", 0.4", 
0.5", 1", 2") and the F606W-band AB limiting magnitude Mini = 
(25, 26, 27, 28), we made 100 images, for a total of 120000. 

In each image, we identify the string-lensed pairs as those 
of the string, and whose separation 



whose centre lies within 87r ^ 



vector is of length < ^T^Gji/c and within 0.05 radians of be- 
ing perpendicular to the straight, simulated string. For the brightest 
magnitude limit, Vnm = 25, this latter condition is relaxed to 0.1 
radians, because the bright, large sources have larger uncertainty 
in their positions. In all we generate 30000 simulated string-lensed 
image pairs. 

Those sources whose centre is further than from the 

String and whose separation is ^ < Sn^^ are chosen as sample 
background pairs, regardless of orientation. We find 500000 such 
background pairs. 

We divide our catalogs of string pairs and background pairs 
along three different axes. First, we split the data into different lim- 
iting magnitude bins (25, 26, 27, 28). The fainter our limiting mag- 
nitude, the smaller the sources we tend to detect. We then further 
divided the pairs into two classes of separation of the string pair, 
"close" and "far," with ^ = 0.5" (a characteristic source size) being 
the cutoff separation. Pairs separated by distances much larger than 
this threshold are well-resolved and deblended, and so should be 
more similar than close pairs. Finally, we divide our sources into 
two classes: "dim" sources that are within one magnitude of the 
limit, and "bright" sources that are more than one magnitude above 
the limiting magnitude. Bright sources tend to have more distinct 
morphology which allows us to distinguish lensed pairs from ran- 
dom pairs more easily. 

We quantify the similarity between sources in pairs using three 
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parameters. The first two are the measured magnitude difference 
between the sources in a pair, 



Figure 1. Example simulated string-lens image pairs caused by a vertical 
string with SttGii/c^ = 1" at redshift 0.5. The ellipses show the ellipticity 
and orientation of each object, as measured with S Extractor. 



Al/= 1^1 - V2\, 
and the pseudo-vector ellipticity dot product, 

E = Si ■ £2 = £l£2 COs(2((/)i — 02)) 
- -Rminor 



(23) 

(24) 
of the 



where si is the S Extractor ellipticity 1 
i^^ source, and is its orientation. 

We also use A0, the angle between the pair separation vec- 
tor and a line perpendicular to the string, to determine how well- 
aligned a string pair is with its string candidate. Note that two well- 
resolved, well- separated and well-measured string-lensed conju- 
gate images will have A(/) — 0, E — 6162, and /SV — 0. 

The distributions of these three parameters are shown in Fig. 2. 
We fit these distributions with the following functions: 

1 

1 



P(A1/) 
P{E) 



AV 



(25) 



Pstring(A(/)) 



Pbg(A0) 



v27r CT 



1 



1 \2 
"2 K/^ct^o) 



27rA(/)o 



2A0n 



(26) 



(27) 



(28) 



We measure separate Vb, ^0 and a for string pairs and background 
pairs; A(/)max is the maximum A(/) we accept (0.05 or 0.1 rad, de- 
pending on Mim). The resulting model distributions are overlaid on 
Fig. 2, and the fit parameters tabulated in Table 1. 

In this table we can see that sources in string pairs have lower 
average /SV and higher E than random pairs. We also see that 
A(/) for string pairs is strongly clustered within | A$| < 0.1. Ran- 
dom pairs have a uniform random A(/) and their distribution is not 
shown. In Fig. 2 we see that our model distributions give good fits 
in all cases. We use these models in the computation of string de- 
tection probability in Section 5.4. 

Incidentally, the term C in equation 27 represents the back- 
ground pairs which are uniformly distributed in < A(/) < 27r. 
We use it to find that there is no more than 5% background pair 
contamination in any of the string pair subsamples. 



5 LENS-FINDING METHODOLOGY 

As explained in Section 2, we look for strings that are straight on 
scales at least as large as 4' (corresponding to > 1.5 Mpc at — 
0.5) by correlating close pairs of faint sources along straight lines 
across small-field images. We give a brief summary of the search 
methodology, before describing it in more detail. 



5.1 Summary 

For each field, we first create a catalog optimised to contain as many 
faint sources as possible, deblending close pairs aggressively. From 
this we produce a catalog of close {0 < 5") pairs of objects. The 
line which bisects each pair is then parametrised by its position and 
orientation in the image — this line represents a putative cosmic 
string. The change of variables between image pair position and 
orientation to bisector line parameters is a modified "Hough trans- 
form:" pairs straddling the same straight line will appear very close 
together in the Hough space (Ballester 1996). We therefore iden- 
tify clusters of points in the Hough space as indicators of possible 
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Figure 2. In each row, we plot the distributions of AV, E and At/) for pairs of faint objects detected in hmiting magnitude = 27 simulated images. 
The rows correspond to the different classes of pairs defined in the text. Top row: close, bright pairs; second row: close, faint pairs; third row: far, bright pairs; 
bottom row: far, faint pairs. 



strings. We then assign a score to each string candidates by evalu- 
ating the similarity of constituent pairs and their alignment with the 
possible string. We visually inspect high-scoring candidate strings, 
remove non-physical image pairs and retest them against a higher 
threshold. In the next two subsections we explain this process in 
more detail. 

5.2 Image pair detection 

We aim to find many sources and not to exclude pairs. We achieve 
both these goals by emulating the catalog production techniques 
in Morganson & Blandford (2008) which are in turn influenced by 
Benitez et al. (2003). Roughly speaking, we used the SEXTRAC- 
TOR program (Bertin & Arnouts 1996) to produce catalogs of faint 



sources. After we make an initial source catalog, we find all pairs 
separated by less than 6 = 5" and make a catalog of pairs. We 
do not probe above 5" because this corresponds to a minimum 
Gfi/c^ = 1 X 10"^ and a typical Gfi/c^ > 5 x 10"^ (includ- 
ing projection and redshift effects) which is not an interesting area 
of parameter space, because it is ruled out with high significance 
by Jeong & Smoot (2007). 

Our source extraction parameters are summarised in Table 2. 
We make three changes to the setup used by Morganson & Bland- 
ford (2008). To make our faint detections very sensitive, we low- 
ered our DETECT.THRESH from 1.7 to 1.6 for WEIGHT type 
weight images. However, we primarily use RMS weight images 
(which were not available for the GOODS data that our cata- 
log techniques were initially designed for), and the two types are 
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Limiting Magnitude 


Subsample 


AVb,str 


AVb,bg 


^0,str 


^0,bg 


Cstr 




A0O 


25 


Close, Bright 


0.350 


0.604 


0.086 


0.006 


0.095 


0.109 


0.046 




Close, Faint 


0.355 


0.511 


0.010 


0.010 


0.126 


0.152 


0.052 




Far, Bright 


0.327 


0.997 


0.057 


0.002 


0.154 


0.157 


0.036 




Far, Faint 


0.382 


0.569 


0.016 


0.000 


0.134 


0.157 


0.044 


26 


Close, Bright 


0.215 


0.707 


0.081 


0.019 


0.131 


0.126 


0.026 




Close, Faint 


0.284 


0.534 


0.041 


0.013 


0.111 


0.153 


0.028 




Far, Bright 


0.285 


0.907 


0.069 


0.004 


0.121 


0.129 


0.021 




Far, Faint 


0.324 


0.633 


0.020 


0.001 


0.125 


0.153 


0.023 


27 


Close, Bright 


0.205 


0.666 


0.086 


0.026 


0.122 


0.145 


0.023 




Close, Faint 


0.234 


0.530 


0.039 


0.014 


0.101 


0.159 


0.027 




Far, Bright 


0.299 


0.942 


0.090 


0.000 


0.125 


0.141 


0.019 




Far, Faint 


0.263 


0.651 


0.033 


0.001 


0.115 


0.150 


0.021 


28 


Close, Bright 


0.200 


0.684 


0.100 


0.023 


0.116 


0.138 


0.022 




Close, Faint 


0.233 


0.545 


0.039 


0.015 


0.093 


0.158 


0.026 




Far, Bright 


0.259 


0.991 


0.113 


0.009 


0.138 


0.144 


0.018 




Far, Faint 


0.249 


0.668 


0.038 


0.004 


0.109 


0.147 


0.020 



Table 1. Fitted parameters for model distributions of the properties of string pairs (str) and background pairs (bg). The close/far dividing separation is 0.5," 
and the bright/faint division is at one magnitude brighter than the limiting magnitude. 



WEIGHT.TYPE 


MAP_RMS 


Filtering FWHM (Pixels) 


1.5 


DETECT.THRESH 


0.895 


DETECTJVLINAREA 


10 


DEBLENDJSTTHRESH 


32 


DEBLENDJVLINCONT 


0.0003 



Table 2. The S Extractor parameters we use to find faint sources and 
close pairs. 



normalised differently. Setting the DETECT_THRESH for RMS 
weight images to 0.895 produces roughly identical catalogs as 
those produced above with weight type WEIGHT. We also changed 
DEBLENDJMINCONT from 0.03 to 0.0003 to promote more ag- 
gressive deblending. 

These changes do not produce many false counts from noise 
peaks in our catalogs. Running SEXTRACTOR on "negative" im- 
ages where the pixels have been multiplied by —1 produces 0.0002 
false sources for every source detected in real images. However, 
complex sources are often deblended into several sources. At this 
stage, the catalogs are inclusive as possible; false pairs will be re- 
jected later in the lens-finding process. 

From our initial faint source catalog, we produce a catalog of 
close pairs. In this catalog, we derive various pair parameters in- 
cluding the centre of each pair, (xo, yo), and the separation vector 
(Ax, A^/). We also compute each pair's AV, E and A$ parame- 
ters (see Section 4.2) and interpolate between limiting magnitudes 
to obtain the statistical properties of string pairs and background 
pairs in Table 1. 

String-lensed image separations can be < ^ < 
SttG/x/c^ sini, so we never set a minimum when look- 
ing for string pairs. However, the number of background 
pairs increases with the separation limit, so we search 
for strings with different maximum 6 to reduce back- 
ground pairs when searching for small separation strings. 
We produce 15 pair subcatalogs for each field with < 
(0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0) 
arcsec. All the work in Section 5.3 and beyond is performed on 
each pair subcatalog separately. 



5.3 The Modified Hough Transform to String Parameter 
Space 

We cannot make a detailed analysis of every possible string orien- 
tation in every field and must focus on straight lines with many 
string lens candidates along them. An efficient way to find sig- 
nificant numbers of aligned pairs is to perform a suitable modi- 
fied Hough transform on all detected pairs, and then accumulate 
a Hough space image (histogram) and search this image for sig- 
nificant peaks. These peaks represent multiple image pairs which 
could have all been lensed by the same straight string. 

To convert a single image pair to its corresponding putative 
lensing string, we find the straight line which bisects the separation 
vector between the objects in the pair. For a pair at mean position 
(xo, yo) and separation vector (Ax, Ay), this line is defined as: 



y = 



Ax 

~^(^ - ^o) -\-yo 

Ay 

[x - XH) + yn 

yn 



{xH^yn) 



(Ax, Ay) 



Ax xq + Ay yo 
Ax2 + Ay2 



(29) 



(30) 



where {xh, yn) are the coordinates of the point where the bisector 
line (string) comes closest to the origin (the centre of the field). We 
convert the impact parameter vector (xH^yn) to polar coordinates 
{Oh, (t>H), because a set of randomly oriented and distributed pairs 
will produce a uniform random field in these coordinates. Oh is 
then the minimum distance (in arcsec) from the bisector line to the 
origin, while (j)H is the polar angle to the point of closest approach 
between the string candidate and the origin. (j)H is measured anti- 
clockwise from the horizontal image axis. We illustrate this trans- 
formation graphically in Fig. 3. This transform is equivalent to a 
standard Hough transform of point sources, except that each pair, 
in addition to having a position in image space, also has an im- 
plied orientation angle. When transformed into Hough space, each 
of these pairs becomes a point source representing a single line 
which goes through a source point rather than an extended Hough 
source representing every line that goes through that source. 

In this new mathematical space, each image pair corresponds 
to a single source. We would like to weight these sources by their 
likelihood of being a string pair rather than a background pair. To 
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Figure 3. We schematically show the transformation from a close pair (in 
physical space) to a potential string (in Hough space). Our pair is repre- 
sented by the two points and S2. The most likely string is the perpen- 
dicular bisector of this pair. Regardless of the location of the pair along this 
string, this string can be defined by the coordinates of the point (xh, Vh) 
where it is closest to the origin (the centre of the data image). When we con- 
vert this impact parameter to polar coordinates (6h,4>h) we can produce a 
Hough space where the points from random pairs are evenly distributed. 




50 100 



Figure 4. Example Hough space image for one of our simulated single 
string-lensed fields. The image extends from — 4"> 9h < 120" along the 
horizontal axis, and — tt < < 3.3 along the vertical axis. The slight 
asymmetries in boundaries allow us to account for various edge effects. The 
simulated string in this field is at (</>, d) = (25^^ 0.0), and can be seen as a 
high- significance object in the the Hough space image (halfway down and 
left of centre). 



this end, we assign each source a non-unit intensity / given by: 



= 1 + 7 log 



Pstring(Ay,^) 

Pbg(Ay,£;) 



\AVo strCrstrA0O 



AV 

sir J AVb bg 
(£;-/Xstr)^ 



(31) 
(32) 



(Note that we approximate P( AF, E) = P{!SV)P{E), since AF 
and E are uncorrelated for all but the brightest, most completely 
copied pairs). By construction, all pairs bisected by the same line 
will appear at the same point in Hough space. If Pstring = Pbg for 
every pair, each pair would have intensity of 1 and the total Hough 
space flux of a feature (string) would be the number of pairs along a 
string. However, pairs of similar sources where Pstring > Pbg are 
weighted more significantly than other pairs. The parameter 7 = 
0.55, and we discuss its purpose and optimisation in Section 5.4. 

After accumulating individual points in Hough space, and as- 
signing them intensities, we convolve the resulting (Bh^ 0//) space 
image by a 2^^ x 4° kernel so that nearby points blend together to 
make multi-pair Hough sources. Four degrees is the typical A(/)i/ 
between the pair separation vector and the actual direction of the 
string. 2" is the typical error in Qh caused by projecting the string 
back towards the origin incorrectly due to A^//. The resulting con- 
volved images contain bright sources corresponding to lines in the 
real space images which have many pairs along them; an exam- 
ple is shown in Fig. 4. We extract bright sources from these im- 
ages using S Extractor with DETECT_MINAREA = 7 and DE- 
TECT.THRESH = 0.16. These settings are optimised to find ob- 
jects with total Hough space flux of 2 units (indicating 2 or more 
pairs of well-matched galaxy images lying along a line in real 
space). Every peak we detect in this Hough space becomes a string 
candidate. 



5.4 Candidate Assessment: Scoring Potential Strings and 
Score Thresholds 

Having eliminated the vast majority of possible string orientations 
and identifying likely string candidates as significant peaks in the 
{Oh,4>h) Hough space image, we examine each candidate more 
carefully using information from each galaxy image pair to assign 
it a score based on its relative likelihood of it being a string lensed 
pair and a random pair. Our first step is to find all pairs which could 
possibly be associated with our string candidate. We cast a wide 
net finding all pairs with (I)h and Oh within 0.1 and 10" of our 
ideal strings. This allows even very incompletely copied pairs to be 
included for consideration as string lensing candidate pairs. 
We then maximise a string "score" defined as: 



S 
Si 



i 

1 + 7 log 

1 + 7 log 
AVi 



+ 



Pstring(AVi, Ej, A(j)i,yi) 

Pbg(Ay,,^„A0,yO 
AVo bg^7bg(/> max bgl/max bg 

TtAVo strCTstr A(/)0 str A^/Q 

AVi , (^z-/ibg)^ 



(33) 

(34) 
(35) 



AVbbg 

{Ei — jjjstr 

2^ 



^ -7log(20, 



2< 



Except for yi and yo, the terms in the above equation are dis- 
cussed in Section 4.2 and 5.3 and represent statistical distributions 
of string pairs and background pairs, yi is the distance of the cen- 
tre of the pair from the string. To simplify the calculation of yi, 
we rotate our coordinates by 00 = Si 0* so that the string angle 
= (/>o + A(/) and each pair angle ^i = + A0i . All the A0's are 
small enough to permit linear approximations. After this rotation, 
Xi is the distance of the pair along the string, is the ^/-intercept of 
the string and yi is the of the centre of each pair from the potential 
string. 

In parameterising our string, we only maximise ^S* with respect 
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to A(/) and 0. All but the last two terms of Eq. 35 are constant, and 
our score is maximised by solving: 



100 



Ef4 



+ 



1 







' A(t) 


















(36) 



At this point, the pair with the most negative Si — correspond- 
ing to a much higher likelihood of being a background pair than a 
lensed pair — is rejected, and and are recalculated until all 
pairs have Si > 0. Only strings with at least three positively scored 
pairs are considered possible strings. This three pair requirement is 
included in our analysis in Section 3. 

Two terms in this algorithm, yo and 7 are not taken directly 
from comparisons of simulated string pairs and background pairs. 
In defining yo we aim to penalize strings which are farther from 
the proposed string than could be allowed by the pair separation, 
its small uncertainty in position or a slight curvature of the string. 
The centroid of legitimate string pairs can be up to 0.5^ (where 6 
is the pair separation) away from the string. The angular size of 
a typical faint galaxy is roughly 0.2" and this leads to a similarly 
sized uncertainty in position when a galaxy image is cut by a string. 
The radius of curvature of a horizon-scale cosmic string should be 
of order 1 Gpc 2.1. 

Across a typical 4' ^ 4 Mpc GO field, a string with a 1 Gpc 
radius of curvature will bend 2 kpc or an observed 0.3" from being 
a straight line. We do not want to penalize strings for being less than 
0.5" or half their pair separation from the string. We therefore set 
yo =max(0.5'' + .56',6'). 

The parameter 7 determines the relative weighting between 
the existence of a pair along a string, and the degree to which the 
pair resembles an ideal string pair. We derive the optimal value of 7 
from our simulations. Our goal is to produce the highest score dif- 
ferential between the string pair and the highest scored background 
pair. We vary 7 in increments of 0.05 and run our algorithm on 15 
simulations with limiting magnitude 27, a string at 2: = 0.5 and at 
string tensions ofSirGfi/c^ = (0.5'', 2") (five each). 7 = 0.55 
produces the highest average S'string - max(5'background). 

Having defined a metric, S, to evaluate string candidates, we 
use simulated string scores to set a threshold score for all potential 
strings. This score will be dependent on the limiting magnitude, 
maximum pair separation and field size (although in practice, the 
field size is fairly constant). We combine these factors into So, the 
expected score for a randomly oriented string in our field: 

^ Ahpsf e ^ c 



(37) 



Where the term ^ooo accounts for the number of resolvable strings 
(the number of PSF areas that fit in our Hough space), and we are 
summing over all the pairs separated by some 6 < ^max which 
has fifteen values for each field as described in Section 5.2. We 
plot the observed scores of all simulated strings with three or more 
pairs versus the expected score in Fig. 5 and determine an empirical 
threshold which accepts 99.5% of simulated strings: 



*S'threshoid(5'o) = max(5.4, 24So'' 



(38) 



Scores will be modified slightly in Section 5.5 after string can- 
didates have been inspected by eye and bad pairs removed, so at this 
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Figure 5. The observed string score versus the expected string score of a 
randomly oriented line for simulated strings. The threshold line is below 
99.5% of observed strings. 



point, we pass all pairs with: 

Si > 0.95'th reshold 



(39) 



We learned by trial that random alignments of two pairs are too 
common to consider as strings, so we also require that each poten- 
tial string have three or more constituent pairs (three pairs detec- 
tions are assumed when we use our results to set limits). 



5.5 Candidate Assessment: Initial Human Inspection 

String candidates which exceed the threshold in equation 39 are 
passed along to more detailed statistical analysis, but first the con- 
stituent pairs of each string are evaluated by eye to exclude false 
pairs from future consideration. The false pairs we exclude come 
in four types: residual cosmic rays, diffraction spikes, large fore- 
ground objects (technically not false pairs, but generally not cos- 
mologically distant) and misprocessed data. We show examples of 
each in Fig. 6. A common problem with these false pairs is that they 
can produce an overdensities along straight lines, thus mimicking 
the cosmic string signal. Strings with false pairs are not automati- 
cally eliminated at this stage unless their total number of valid pairs 
goes below three. 



5.6 Candidate Assessment: Detailed String Analysis 

Having eliminated bad pairs from the potential stings, we can ex- 
ecute two more rigorous cuts knowing that our results will not be 
significantly influenced by false pairs. The first cut is to require 
that the string go between each pair (or nearly so) while allowing 
for some string curvature. The second cut is to require that the score 
per pair is similar to that observed in simulation. 

To fit a curved string to a set of pairs, we first fit a straight 
line using Eq. 36. We then rotate the pairs by —Acj) so that they 
are perpendicular to this fitted line. We define a Cartesian system 
in which x is along the linear approximation of a string candidate 
and y is perpendicular to that line. We fit the pair positions with: 



■ 6 + axi 



(40) 



We intentionally do not fit the curve with three degrees of free- 
dom at once. The spatial terms Xi are centered so that Xi = Ois in 
the middle of the string. Fitting the quadratic term in equation 40 
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Figure 6. The four basic types of false pairs rejected during human inspection: cosmic rays not removed by our automated algorithm (upper left), a diffraction 
spike (upper right), a foreground cluster (lower left) and a hne of bad pixels due to poor field processing (lower right). 



allows for symmetric curvature across the field, as we would expect 
for a string with a very large radius of curvature. A full three-term 
quadratic fit has the undesirable property of allowing the fit to curve 
out sharply at either end. 

As noted in Section 2.1, the curvature of the string is believed 
to be small. We quantify this by assigning the following Gaussian 
prior on a: 

P{a) oc e (41) 



The most probable curved string given the data is then given by: 




Once this fit is obtained, we examine all pairs in the field that this 
curve does not pass through and whose closest source is more than 
0.2" away from the curve. The pair whose closest source is farthest 
from the curve is excluded. The fit is recalculated, first as a line and 
then as a curve until all pairs are either split by the curve or within 
0.2" of being split by the curve. This more strict filtering elimi- 
nates many high separation string candidates, but does not affect a 
significant number of simulated straight strings. 



Our final method of string candidate rejection is to assert a 
cutoff on string score as a function of the number of pairs associated 
with the string, S{n). In Fig. 7, we plot S(ri) for our simulated 
strings and find that 

OThreshold(n) = — (43) 

1 + n 

is a threshold that accepts more than 99.5% of strings with 3 or 
more observed pairs. 

String candidates that survive these two cuts are passed along 
to a final round of field-dependent inspection in which neighboring 
fields in surveys and multi-band images are used (where possible) 
to eliminate string candidates. These methods require more human 
involvement and are not generally applicable, so we present them 
separately in Section 7. Having outlined the end-to-end string can- 
didate generation procedure, we next describe the imaging dataset 
we use in our string search. 



6 SURVEY DATA: THE HST/ACS ARCHIVE 

We now describe the input data for our cosmic string search. As 
previously discussed, we use HST optical imaging data for its high 
resolution and background source density. For ease and homogene- 
ity of processing we focus on images taken with the Advanced 
Camera for Surveys(ACS). The total sky area imaged with ACS is 
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Figure 7. The observed string score versus the number of pairs. The thresh- 
old hne is below 99.5% or observed strings. 
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Figure 8. HAGGLeS HST/ACS field positions, in galactic coordinates. The 
symbols are not to scale. 



Survey 


No. of fields 


Area 


Mean depth 






(deg2) 


(AB mag) 


GO-S 


318" 


1.08 


25.8 


GO-M 


346" 


1.18 


26.3 


COSMOS 


575 


1.84^ 


26.5 


AEGIS 


63 


0.20 


25.8 


GEMS 


63^ 


0.18 


25.3 


GOODS 


35 


0.10 


28.6 


Overall: 


1400 


4.48 


26.2 



Table 3. HST/ACS archive data used. See text for our source detection pa- 
rameters, and the definition of "depth." Notes: a) This is the number of 
fields remaining after visual inspection for high cosmic ray or star density 
(Section 6.2). h) A pointing overlap area of 4% was assumed for each of the 
large surveys, as in Scoville et al. (2007). c) For the purposes of this paper, 
GEMS does not include GOODS-S area: we have subtracted the overlap 
area (which we estimate to be about 0.015 deg"-^, Rix et al. 2004) between 
the two surveys from the GEMS area. 

approximately 12 square degrees. However, we require the images 
to be sufficiently deep and high galactic latitude to be able to detect 
significant numbers of faint galaxies at redshift 1-2, and for those 
images to be largely free from confusing cosmic rays. We selected 
4.5 square degrees of ACS imaging data that met these criteria (al- 
though see Section 6.2 below for further discussion of cosmic ray 
contamination); about half of this data comes from the large pro- 
grams GEMS (Rix et al. 2004), GOODS (Giavalisco et al. 2004), 
AEGIS (Davis et al. 2007) and COSMOS (Scoville et al. 2007). 
The remainder are data from extragalactic Guest Observer (GO) 
observations. For reasons explained in Section 7 we divide the GO 
images into GO-M (fields with imaging in multiple filters) and GO- 
S (single filter fields). Approximately 0.75 square degrees of the 
GO-M imaging had more than 2000 seconds exposure time, and 
was made public before November 2005: this subset was searched 
for galaxy-scale strong gravitational lenses by Marshall et al. (in 
preparation) as part of the HST Archive Galaxy- scale Gravitational 
Lens Search (HAGGLeS). Here, we use a larger set that includes 
more recently-observed fields, fields that were only observed in one 
filter and fields with a little as 1500 seconds exposure time. We re- 
tain the need for 3 exposures per field for effective cosmic ray re- 
jection. We summarise the HST/ACS imaging data used in Table 3. 

The geometry of the GO survey is shown in Fig. 8. The avoid- 
ance of the galactic plane area can be seen, but otherwise the 



HST/ACS archive does a reasonable job of uniformly sampling the 
sky. 



6.1 Survey field depths and areas 

We initially search through individual fields from large surveys as 
though they were independent entities. When determining the lim- 
its implied by null detections in each survey, however, we view 
each contiguous survey as one large area (with GOODS North and 
GOODS South being two areas). This method automatically ac- 
counts for field overlap in surveys, reducing areas by 4% (as ap- 
propriate for the largest survey, COSMOS, Scoville et al. 2007). 
We could search for strings across full surveys as though they were 
individual fields, but our measurement of Oh becomes less accu- 
rate as fields get larger. In Section 7.1, we discuss how we use the 
contiguous nature of large surveys to eliminate string candidates. 
The GO fields are predominantly single ACS pointings, although 
in some cases we mosaiced overlapping fields together. 

Very few instances of GO field overlap remain. We account 
for overlaps between the large program fields by working with a 
single catalog of unique sources for each large program area; we 
then correct the large program survey areas downwards by 4% (as 
appropriate for the largest survey, COSMOS Scoville et al. 2007) 
when computing string densities and parameter constraints. 

The "depth" of each field is a label allowing the sensitivities 
to source detection to be quantified in each field. It is defined as 
the F814W AB 5-sigma limiting magnitude for an extended source 
measured within a 0.4" radius circle, were it a faint blue galaxy at 
redshift 1.5. This definition was adopted to allow comparisons be- 
tween fields observed in different filters: specifically we assume the 
Scd template spectral energy distribution of (Coleman et al. 1980) 
for the transformations. The depths we derive match well the ex- 
tended source F814W AB limiting magnitudes calculated by the 
AEGIS project team (Davis et al. 2007) and the 90% faint galaxy 
completeness limit estimated for their images by the COSMOS 
team (Leauthaud et al. 2007). We chose the F814W filter for con- 
sistency with the model redshift distribution of (Leauthaud et al. 
2007) that we have adopted. This is plotted for each of the surveys 
in Fig. 9. 
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Figure 9. Source redshift distributions assumed for the different surveys. 
The "GO" curve shows the mean redshift distribution: in practice, we used 
a different distribution for each field according to its depth. 
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Figure 10. The distribution of brightest pixel intensity distribution in field 
ACSJ0024 18-020750. We see a characteristic double Gaussian distribution 
with the brighter of the two peaks being due to cosmic rays. 



6.2 Image processing and Cosmic Ray rejection 

In the case of the large programs, we used the high level archive sci- 
ence product images provided by the project teams. The GO images 
were reprocessed from the uncalibrated data using the multidrizzle 
pipeline developed for the HAGGLeS project (Marshall et al. in 
preparation). This processing included a visual inspection step to 
remove cosmic ray clusters, satellite trails, scattered light and other 
confusing artifacts. 

Many of the GO fields were produced from only a few raw 
images. While HAGGLeS aims to produce uniformly high-quality 
data, many fields suffer from cosmic ray corruption not recorded 
in the weight image. Since we are not interested in perfectly accu- 
rate catalogs but rather rare, high-level correlations, we use a sim- 
ple method to extract the vast majority of these sources and accept 
some impurity. 

We have many processed images with false sources caused 
by cosmic rays. Typically, these detections are around the edges of 
images where the input raw images did not overlap, but in some 
fields, cosmic rays corrupt the image interior as well. The HAG- 
GLeS weight images usually mark these pixels as having high vari- 
ance, but this process is likely to fail when there are fewer input ex- 
posures. Fortunately, cosmic ray events tend to deposit more energy 
on a pixel than all but the bright sources. Examining the intensity of 
the brightest pixel in every detect source in a typical HAGGLeS GO 
field produces a double Gaussian distribution as shown in Fig. 10. 

The higher intensity peak is due to cosmic rays, while the 
lower intensity peak is due to astronomical sources. We fit the dis- 
tribution with two Gaussians centered at log / = /ii and /i2 and 
with widths ai and (72, and reject sources with log / > /io where: 

/il + 2(71 + /i2 - 2(72 r ^ o o .... 

Mo = , for fii < fi2 - 2(71 - 2(72 (44) 

Ml + CTi + /i2 - (72 r 

fio = , for 

M2 — 2(7i — 2(72 < Ml < M2 — CTl — 

We do not filter out bright sources if the peaks cannot be resolved 
at the (71 + (72 level. Generally, these fields have very low cos- 
mic rays counts, but some fields are dominated by cosmic rays and 
were excluded during visual inspection. After we have identified 
a population of cosmic rays, we search for significant overdensi- 
ties around the edges of our image and remove these edges if an 
overdensity exists. This filtering is imperfect since the cosmic ray 



intensity peak has a long negative tail and the source peak inten- 
sity has a long positive tail, but we estimate roughly 95% purity 
in typical samples. We eliminated 74 of the 738 GO fields as being 
unsuitable for string detection due to excessive cosmic rays, stars or 
foreground clusters (which prohibit string detection just like excess 
cosmic rays) in this way. 



7 RESULTS 

The search methods described in Section 5 are successful at greatly 
reducing the number of potential strings in our surveys, but we must 
address the small number of remaining candidates in a rigorous way 
in order to produce detections or robust upper limits on string ten- 
sion and concentration. Our methods sort through millions of pos- 
sible strings in 1400 fields, leaving only 50 string candidates which 
we must study in more depth. These remaining candidates show ev- 
ery indication of being the statistical tail of coincident alignments. 

In Table 4, we show how many string candidates existed at 
each stage of our analysis. Each of our fifteen Hough images can 
resolve roughly 1000 strings, so we are in effect examining 15000 
potential strings per field. Only 1638 potential strings (1.2 per field) 
pass initial automated inspection described in Section 5.4. These 
candidates were all examined by eye and had any obvious bad pairs 
removed. Only 50 of the remaining candidates passed the final au- 
tomated test in Section 5.6 with at least 3 pairs and a score per 
pair consistent with simulated strings. The 3 -pair requirement is 
particularly effective at eliminating potential strings in bright limit- 
ing magnitude fields (AB Mag < 26) and small separation strings, 
because there are fewer background pairs in these cases. Poten- 
tial strings in high limiting magnitude fields (AB Mag > 27), and 
those that produce high angular separations, have higher thresholds 
which require more chance alignments, and so are also efficiently 
excluded. It is the fields with intermediate limiting magnitudes (26 
< AB Mag < 27) that produce the most automatically-assessed 
candidates, and we have relatively few candidates remaining from 
the GO survey (which has many shallow and deep fields) and none 
from GOODS (which is uniformly deep). 

The numbers of pairs in the 50 remaining string candidates 
suggest that they are false detections. While we probe separations 
in the range 0.3"< < 5", we only find candidates with maxi- 
mum separations of 0.6" or greater. If a string were making 0.6" 
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Survey 


Number of Fields 


Resolvable Strings 


Passed Assessment 1 


Passed Assessment 2 


Possible Strings 


GO-S" 


318 


4.6 X 10^ 


323 


4 


4 


GO-M^ 


346 


5.3 X 10^ 


382 


8 





COSMOS 


575 


8.7 X 10^ 


631 


31 





GEMS 


63 


1.0 X 10^ 


135 


5 





AEGIS 


63 


1.0 X 10^ 


120 


2 





GOODS 


35^ 


0.3 X 10^ 


46 








Total 


1400 


21.9 X 10^ 


1637 


50 


4 



Table 4. The number of fields from each survey of HST data, the approximate number of strings we can resolve in our Hough space, the number of candidates 
which pass the assessment described in Section 5.4, the number of candidates which pass the assessment described in Section 5.6 and the number of candidates 
that we cannot reject. Notes: a) The GO-S fields are the GO fields with only a single color filter of data, h) The GO-M fields are the GO fields with multiple 
color filters taken, c) The number of GOODS fields is the number of panels produced by the GOODS team and does not represent individual HST exposures. 



separated image pairs of high redshift sources, we would expect it 
to make roughly: 

0.5 X 5000 (9max 



< rip 



= 7.2 



(45) 



210^^ 

Where the factor 0.5 is to account for the fraction of sources which 
are significantly behind the string, 210" and 5000 are typical ACS 
width and the number of sources in a panel of our shallow sur- 
veys, respectively. The majority of our string candidates exhibit 
larger separations, indicating that we would expect even more than 
7 pairs. But of our 50 string candidates, 38 consist of 3 pairs, 11 
consist of 4 pairs and 1 consists of 5 pairs. 

The above argument is not quantitative enough to explicitly 
figure into our analysis. But we can use other objective methods to 
show that string candidates in large surveys and fields imaged in 
multiple filters do not represent real cosmic strings. 

7.1 Detailed string candidate inspection 

The 50 candidate strings described above were selected by geome- 
try and single filter image pair morphology alone. To test the string 
lens hypothesis further in each case we use different methods de- 
pending on whether the string candidate in question is in a large 
survey or an isolated field, and whether it has been imaged with 
multiple filters or just one. After all string candidates are exam- 
ined, we can eliminate all but four candidates using objective cri- 
teria, and we believe that even these final four candidates could be 
eliminated with more data as we discuss in Section 7.2. 

In large contiguous surveys, our assumption of long, straight 
strings requires that any strings detected in a single field be de- 
tected as collinear strings in neighbouring fields. We illustrate this 
idea for COSMOS in Fig. 1 1 by plotting the extrapolation of every 
string candidate into the rest of the survey. While there are some co- 
incidentally aligned strings, there is no pattern of strongly collinear 
string candidates to indicate a real cosmic string. To make this anal- 
ysis more precise, we use the techniques in Section 5.3 to make a 
single survey-wide Hough image of every string candidate in the 
survey (Fig. 12) and see if there are any significant alignments. Our 
measurements of string candidate orientation angle are precise to 
0.01 radians in a field, but we want to allow for modest curvature 
so we look for pairs whose (j)H are within 0.1 of each other and 
whose Oh are within 0.14° (allowing the uncertainty in (j)H to be 
projected across the 1.4° survey). Only one pair of string candi- 
dates survives this cut. This pair is composed of two candidates in 
acs_100036+0205_unrot, which are actually parallel but separated 
by more than an arcminute, rather than being collinear. Analogous 
Hough space images for GEMS and AEGIS also show no align- 
ment of string candidates. 

This argument breaks down for any string candidates that cut 




150.8 150.6 



150.4 150.2 150 149.8 
RIGHT ASCENSION 



149.6 



Figure 11. The 31 string candidates in COSMOS with the location of the 
string at the foot of the vector and the string projection along the line of 
the vector. The graphing area is roughly the same as the COSMOS survey. 
Every string candidate, when extrapolated, would cut across a significant 
fraction of the survey and so should thus be detected in multiple fields. 
There is no strong pattern of string candidates in a line to indicate this. 



across the corner of a large survey, or any that cut across a fairly thin 
survey. Fortunately, none of our string candidates cut across single 
corner fields. The survey string candidate least likely to be observed 
in neighboring fields is the second string detected in the EOS- 16- 
02. The AEGIS survey is 3 fields wide and 21 fields long, and the 
projection of this string cuts across only three fields. This string 
candidate consists of three 1.0" pairs and the typical AEGIS field 
contains 5000 sources, so, analogous to Eq. 45 it should produce 
an average of 11.9 events per field as above. The chance of such a 
string producing only three events in one field and less than three 
events in two neighboring fields is essentially zero. 

Eliminating string candidates across the GO fields requires a 
different approach, because the GO fields are not contiguous. For 
the GO-M fields, we can eliminate candidates by finding pairs for 
which one member of the pair has significantly different colour than 
the other. String lensing is ideally an achromatic effect, but it is 
possible that a galaxy could be incompletely copied so that, per- 
haps, the blue star-forming core only appeared in one image, and 
the two images could have very different colors. For this reason, 
we do not automate the process and search by eye for sources with 
different colour that did not seem to be cut in half. We eliminated 
enough pairs from each potential string to bring the pair total below 
three. We show an example of a colour-rejected pair in Fig. 13. This 
multi-band analysis could reasonably be automated, and points to 
the usefulness of multi-band analysis in future string searches. 
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Figure 12. The Hough image of potential COSMOS strings. We only find 
one near overlap of string candidates. Examination of the candidates in de- 
tail shows that they are parallel lines detected in the same field and separated 
by more than an arcminute. 




Figure 13. Example multi-band dropout in the GO field ACSJ065819- 
555630. The F814W-band detection image shows a pair of two sources 
(top), but in the F435W filter (bottom) only one source is apparent. 



7.2 The four remaining string candidates 

We are left with four potential strings from the single band GO 
fields (see Table 5 for more information). We emphasise that there 
is nothing particularly suggestive about these remaining candidates, 
but only a lack of data that prevents us from eliminating them. We 
display the constituent pairs of each string candidate in Fig 14. 
These string candidates come from some of our noisiest data, and 
the images we display here are scaled differently from other data. 
Even qualitative analysis of these images suggests that these pairs 
are not exact copies. We show in this section that statistical analysis 
strongly suggests that these string candidates are indeed false de- 
tections. But we also show the limits on Gfi/c^ that each of these 



string candidates would imply if it were the result of a real cosmic 
string. 

The rejection of the 46 potential string candidates which met 
the same criteria as the remaining four allows us to estimate the 
probability that these remaining candidates are false detections. If 
we start with the prior assumption that the fraction of string can- 
didates which are false detections, /, is a random variable cho- 
sen uniformly between and 1 (a generous assumption) and label 
the condition of having 46 false detections "46 false", then Bayes' 
Theorem gives us a distribution on f: 



P(/|46 false) 



P(/)P(46false|f) 
P(46 false) 



(46) 



And in turn, the probability of the remaining three detections being 
false is: 



P(4 false) 



/' 

Jo 



47f'*''f'^df' = — = 0.92 
51 



(47) 



This 0.92 is a minimum probability that makes use of a gen- 
erous prior and does not take note of the fact that the scores and 
numbers of pairs found in the four remaining string candidates are 
atypically small for cosmic string detections. We thus believe that 
these detections are spurious, although we cannot prove it with cur- 
rent data. 

If we accept these detections as genuine, we can calculate the 
implied tension and density of cosmic strings as shown in Table 5. 
We show the limits (as calculated in Section 3.3 for potential strings 
in Fig. 15. 

From Table 5, these four strings would imply a most likely 
G/i/c^ of 6.5 X 10"'^, 1.0 X 10"^ 1.5 x 10"^ and 2.4 x 10"^ 
respectively and a most likely Qs of 7.7 x 10~^, 1.1 x 10~^, 
1.5 X 10~^ and 2.2 x 10~^, respectively. Any of these strings 
would be consistent with limits from other direct detection meth- 
ods, but inconsistent with the stricter indirect CMB detection limits 
(see Section 3). 

Following from the results of the previous section, it would be 
useful to image these four string candidates with comparably high 
resolution imaging in a second filter. In the absence of such data, 
we will present our limits on string concentration and tension in the 
next section with the multi-band GO (for which we have no string 
candidates) and total GO as separate limits. 



8 A NEW DIRECT DETECTION UPPER LIMIT ON THE 
STRING TENSION 

In this section we use the results from our direct detection pro- 
gramme to set direct detection limits on the cosmic string tension, 
Gii/c?, and string density, Qs, which vastly improve upon previ- 
ous optical direct detection limits and are competitive with indirect 
detection methods such as CMB power spectrum analysis. This im- 
provement is largely due to our using the many scattered GO fields 
rather than just one or two contiguous surveys. Our limits become 
much tighter if we assume that the four detections in the GO single 
band images are false detections. 

In Figure 16, we show the contours in Gii/(? — ^ space for 
which we would detect a string with 95% probability, derived from 
each of a range of HST/ACS surveys of increasing scale. We use 



16 Morganson et al. 



Field 


^max 


Right Ascension 


Declination 


Slope 


N ■ 

^ ^ pairs 


Threshold 


Score 


logl0(G/i/c2) 


logio(^s) 


ACSJ195021-405350 


1.0 


297.59409 


-40.88897 


0.1717 


4 


9.43 


10.25 


-6.19 ±0.27 


-5.11 ±0.29 


ACSJ181424+411240 


1.4 


273.60678 


41.20152 


0.6494 


3 


6.58 


7.88 


-6.00 ±0.26 


-4.96 ±0.28 


ACSJ0421 15+193600 


2.0 


65.31551 


19.59700 


0.1949 


3 


7.09 


9.11 


-5.81 ±0.21 


-4.81 ±0.23 


ACSJ144146-095150 


3.0 


220.44917 


-9.86378 


0.8195 


3 


5.55 


5.70 


-5.628 ±0.051 


-4.65 ±0.10 



Table 5. The four candidates which passed two rounds of automated analysis, were not part of large surveys and lacked multifilter data. These candidates are 
all in the single filter GO survey (GO-S). 





Figure 14. Example string-lens image pairs from field ACSJ195021-405350 (left), ACSJ181424+411240 (left-centre), ACS J042 115+193600 (right-centre) 
and ACSJ 144 146-095 150 (right). Each image is rotated such that the most likely string is vertical. The ellipses show the ellipticity and orientation of each 
object, as measured with S Extractor. The quality of these images was, on average, fairly poor, and they had to be scaled differently from other data images 
so that sources we visible above the noise. This accounts for the different appearance. 



the formulae from Section 3. We also show the contours for our 
combined dataset, the projected contour if we doubled its size, and 
the contour that would be obtained by a putative Joint Dark Energy 
Mission (JDEM) weak lensing survey covering 10000 square to an 
AB limiting magnitude of 26 (Aldering et al. 2004). These contours 
are not confidence intervals and make no prior assumption on the 
distribution ofGii/c? or ^. 

Figure 16 illustrates the sharp dependence of our ability to 
probe G^/(? on the number of fields being searched through. This 
dependence is due to the requirement that a moderate redshift string 
cross through at least one field. At the theoretically-predicted val- 



ues of ^ = 60, strings at observable redshifts would cross fewer 
than 1% of the ACS fields. Very high Gji/c? strings would be ob- 
served nearly any time they cross a field, but moderate redshift 
strings with smaller Gfi/c? become undetectable if their projection 
along the line of sight, sini, is small. When we add the 318 single 
filter GO fields to our analysis, we are exponentially decreasing the 
chance that every cosmic string "missed" every one of our fields. 
This leads to a marked decrease in our limits on Gfi/c^ . 

We produce proper confidence limits using the methods in 
Section 3.3. As a reminder, the prior onGfi/c? is log-uniform from 
10~^ to 10~^, and the prior on ^ is Gaussian in logspace center- 
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Figure 15. The G^jc^-^ limits if we assume that our string detections in ACSJ195021-405350 (top left), ACS Jl 8 1424+4 11240 (top right). 
ACS J042 115+193600 (bottom left) or ACSJ 144 146-095 150 (bottom right) is genuine. The inner contour is the 68% confidence limits and the outer con- 
tour is the 95% confidence limit. Our priors are described in Section 3.3. 
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Figure 16. The G^xjc^-^ 95for GOODS, AEGIS, GEMS, COSMOS, GO multifilter and GO total surveys (left). We show these contours for all sur- 
vey s+multiband GO, all survey s+GO, doubhng the data from all survey s+GO and a proposed IDEM weak lensing survey (right). 



erd around log(^) = log(60) and with width log(2)/2. The limits 
for each GO field were calculated independently, while the limits 
for each survey (with GOODS North and South as separate sur- 
veys) were calculated as single large fields. Limits for individual 
GO fields and surveys were then combined. 



In Figure 17, we show the 68% and 95% confidence contours 
in string G/i/c^ — ^ space using all surveys and the multifilter 
GO fields. We also show the analogous contours calculated using 
surveys and all GO fields. For surveys and multifilter data, our 95% 
confidence limits are G/x/c^ < 2.3 x 10"^, logio ^ = 1.72±0.41 
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Figure 17. The 68% and 95% confidence limits for all survey s+GO-M (left) and all survey s+GO (right). Our priors are described in Section 3.3. 



and Qs < 2.1 X 10 ^ . Using all of our data, our 95% confidence 
limits are G/i/c^ < 6.5 x 10"'^, logio ^ = l-^O ± 0.30 and 
Q.S < 7.3 X 10~^. Our limits on ^ are particularly sensitive on the 
prior we take in Eq. 22, but we have lowered the expected value of 
^ from 60 (the prior) to 52 (using GO-M) or 32 (using all data). 



9 DISCUSSION 

Our string constraints could be lowered somewhat by new data that 
will soon be delivered by several different projects. The recently de- 
ployed HST Wide Field Camera 3 will provide 7 square arcminute 
exposures of high resolution data with similar limiting magnitudes 
to the Advanced Camera for Surveys (Kimble et al. 2006): this new 
GO data, combined with the growing ACS imaging archive should 
provide at least a factor of two increase in sky coverage. Looking 
ahead, the James Webb Space Telescope (JWST) is scheduled to 
launch in 2014, and will provide more images of similar depth and 
resolution to HST in the infrared (Gardner et al. 2006), albeit over 
smaller (:^ 4 square arcmin) fields of view (and which may already 
have been observed in the optical with HST). These data would 
certainly enable the constraint projected in the "doubled dataset" 
curve in the right-hand panel of Figure 16. 

Can we envision extending this study to much larger surveys? 
As we showed in Paper I, large ground-based projects like the Dark 
Energy Survey (The Dark Energy Survey Collaboration 2005), the 
Panoramic Survey Telescope and Rapid Response System (Pan- 
STARRS) (Jedicke et al. 2007) and the Large Synoptic Survey Tele- 
scope (LSST)(Ivezic et al. 2008) would provide sufficient area and 
source counts for string searches, but lack the angular resolution 
to further constrain Gfi/c^ . Weak lensing surveys have been pro- 
posed for the Joint Dark Energy Mission (JDEM) that cover 10000 
square degrees to an AB limiting magnitude of 26 (Aldering et al. 
2004), although the exact details of the JDEM mission are still in 
flux (Gehrels & Albrecht 2009). 

Searching through more data will indeed lower our detection 
limits. If we were to double the number of HST fields we searched 
through, our 95% confidence limits would fall by a factor of 2-3 to 
G^/c^ < 2.5 X 10"^ and Qs < 3.1 x 10"^ The proposed JDEM 
survey would lower these confidence limits by a further factor of 
just 1.7 to Gfi/c^ < 1.5 x 10~^ and Qs < 1-8 x 10~ 6: a marginal 
improvement over the doubled HST limits, given the 3 orders of 



magnitude increase in survey area! At this point, resolution, and 
not lack of data would be the limiting factor. In fact, we argue that 
with this work we have within sight of the realm of diminishing 
returns, with only the next factor of two increase in survey scale 
gaining us a factor of two in parameter constraints. 

The computational and human requirements for this work are 
modest and could be easily scaled up to analyse the forthcoming 
increased HST archive dataset. We did not use any special compu- 
tational tools other than parallel batch processing, and all compu- 
tation times quoted below are for modest desktop processors. We 
produced 6000 full HST-size simulated string images and 120000 
10" wide string images (see Section 4) automatically in one month 
of processing time. The first round of automatic analysis (see Sec- 
tion 5.4) was performed on simulated data in one week and on all 
HST data in three days. The manual filtering described in Sec- 
tion 5.5 was performed on automatically-generated webpages for 
each field and required an average of twenty seconds per poten- 
tial string for a total of ten hours of human work. In the future, 
the extraction of cosmic rays and other false sources would need 
to be automated for JDEM-scale datasets; the greater uniformity 
of such a survey's images would make this process easier. The fi- 
nal automated phase described in Section 5.5 took essentially no 
computer time. The analysis of the remaining candidates (see Sec- 
tion 7.1) took several hours and would have to be automated for 
larger projects. 

The technique for long string detection developed here will 
continue to be applicable to wider searches for lower tension, 
smaller image separation, and more sparsely distributed strings. 
Previous direct detection methods relied on finding single distinc- 
tive string lensing events (Sazhin et al. 2007) or searched for a 
large excess of similar pairs (Christiansen et al. 2008). The for- 
mer method requires bright sources at large separation and thus 
does not make full use of survey depth or resolution. This prevents 
it from probing the smallest Gfi/c^ . In addition, it is difficult to 
characterise a single pair as being caused by a string lensing event. 
The latter method does not efficiently extract string lensing events 
from the background and is only effective if the field is crowded 
with strings. By searching for several pairs along straight lines, we 
probe to smaller separations (which imply smaller Gji/c^), exam- 
ine fainter sources (whose exponentially large number allow us to 
probe smaller ^) and limit the possibility that a real string will be 
obscured by background pairs. 
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However, there are several straightforward improvements one 
could make to our string searching methods that could improve 
both the final limits and the searching efficiency in future projects. 
Firstly, our pair characterisation, based on magnitude, ellipticity 
and the alignment of a pair with the proposed string, does not use 
all the information available, and is unable to recognise a very in- 
completely copied pair. One could instead match the images of 
each pair pixel by pixel (as performed by Agol et al. 2006), and 
allow incomplete copying to account for galaxy images "cut" by 
the string. We did not not implement this method on the grounds 
that our search focuses on barely resolved pairs for which there is 
little information, and because it would vastly increase computa- 
tional time. One could imagine a hybrid survey whereby candidate 
image pairs were fed to an automated pixel-matching routine; this 
would both improve the purity of the candidate samples and also 
reduce the human inspection time of the search. 

Using exclusively multifilter data would also make searching 
more efficient. None of our automatic string detection algorithms 
use color information because both COSMOS and half of the HAG- 
GLeS GO fields are single filter. But we saw in Section 7 that 
checking for color consistency across pairs effectively eliminates 
potential string candidates. It may be impossible to guarantee mul- 
tifilter data for large, diverse archival datasets like the GO fields, but 
the kind of homogeneous survey like IDEM would automatically 
provide color information which could be profitably exploited. 

The main reason why our limits on string concentration and 
tension are much stricter than those of previous works is our use 
of a survey comprised of many scattered small fields. Because the 
probability of a string crossing a field scales as the linear size of the 
field and not its area, searching for cosmic strings in randomly scat- 
tered fields is more effective than searching many contiguous fields 
in a survey. This method also has the practical advantage that it 
does not require a dedicated survey and naturally thrives on archival 
data. Essentially any high resolution extragalactic images without 
overwhelmingly large objects in them can contribute to our lim- 
its on cosmic strings. However, the relatively small size of the ACS 
field means that at the low lensing cross-sections of interest we only 
expect to see a small number of lensed image pairs per field. This 
may be a problem for string searches based on the smaller WFC3 
and JWST fields of view. 



10 CONCLUSIONS 

From an exhaustive search of some 4.5 square degrees of archival 
HST/ACS imaging data, we have derived significant direct detec- 
tion limits on the tension and density of cosmic strings, employing 
a variety of techniques that could be applied to future searches. We 
draw the following conclusions: 

• Using multi-filter data from the HAGGLeS, COSMOS, 
GEMS, AEGIS and GOODS surveys, we are able to constrain the 
dimensionless string tension Gfi/c? to be below 2.3 x 10 ~^ and 
the energy density in long strings to be below 2.1 x 10 ~^ with 
95% confidence. 

• Extending the search to single filter GO imaging, we found 4 
string candidates that we were unable to reject formally on the basis 
of the data in hand. Applying our experience with the multi-band 
data, we suggest that these detections are false positives: under this 
assumption we find stronger upper limits oiGji/c? < 6.5 x 10~^ 
andQ, < 7.3 x 10"^ 

• Unlike previous optical imaging searches, we cover sufficient 
sky area to make these limits universal: if there were cosmic strings 



present in the redshift < 1 universe with tension and density greater 
than our limits, we would have detected them. 

• Color information was found to be important in rejecting 
string candidates: high resolution imaging in multiple filters would 
be the most efficient way of ruling out the candidates presented 
here, and any future single-filter detections. 

• Using the technology developed in this study, the upcoming 
factor of two increase in sky area imaged at comparably high reso- 
lution expected from the refurbished HST should enable these up- 
per limits on string tension and density to be further reduced, by 
slightly more than a factor of two. 

This work represents the the first direct detection limits on 
cosmic strings that are both competitive with and complimentary 
to other indirect methods. Our limit on the string tension is only 
a factor of 2-3 higher than the indirect upper limit of Gfi/c^ = 
3 x 10 ~^ derived from the CMB power spectrum. We have shown 
that direct detection of strings by their gravitational lensing effect 
has the advantage of constraining both G/i/c^ and ^, and does not 
require that strings emit gravitational radiation, as the pulsar timing 
constraints require. Perhaps most importantly, direct detection is 
the only method currently available that would provide the precise 
location of a cosmic string for future study. Prospects for using near 
future data to improve upon these limits are bright; however, we 
anticipate that this will bring us into the era where large increases in 
survey area would bring relatively small changes in the limits. The 
next significant advance in this field maybe the increase in angular 
resolution promised by future radio surveys. 
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APPENDIX: THE PROBABILITY OF A 1 RADIAN 
STRING CROSSING A FIELD 

When calculating the probability of a rectangular, Oi x O2 field ly- 
ing on a projected 1 radian string such that Oc string crosses the 
field, we start by making useful approximations and other simpli- 
fications. We use the fact that our fields are much smaller than a 
radian so that the sky can be approximated as a rectangle with area 
47r steradians. We ignore the end effects and curvature of the string. 
As shown in Fig. 18, we imagine that the centre of our field is a dis- 
tance from our string, and our field is oriented with the angle 
where = implies that the long sides of our field (with length 



20 Morganson et al. 



If a string crosses along the long side of the field: 




Figure 18. A string crossing a x 62 field. The distance from the centre 
of the field to the string and orientation of the field are 6 and 0, respec- 
tively. The length of string crossing the field is Oc. The lowest corner is 
Co. The adjacent corners separated from Co by 61 and 62 are Ci and C2 
respectively. 



O2) are perpendicular to the string. We need only consider the cases 
where ^ > and < < 7r/2. The rest of parameter space is 
proportional by symmetry. The probability of a string crossing our 
field with toal overlap Oc is then: 

Prad(^c) ^ \ r dO # 5{e',{e, 0) - Oc) (48) 

Jo Jo 

From Fig. 18, we can see that the vertical displacement be- 
tween each of the three corners and the string is: 



Co 



Cl 



C2 



■ cos((/) - 0o) 
■cos(0 + 0o) 
^ cos((/) + 0o) 



(49) 
(50) 
(51) 



Where Oo = y^Oj + 0^ and 00 = arctan(6>i/6>2). 

A string can cross our field in three different ways: along the 
short axis (^1), along the long axis (O2) or cutting a corner. If the 
string crosses our field along the short axis (entering and exiting 
via the two long sides), it's crossing length is: 



i— -, for < 7r/2 - 00 and 6 < a (52) 



cos(0) 



Using Eq. 48 we can convert this into the component of Prad that 
represents the string crossing the field along the short axis, Pshort : 



=) 



for 6>i < 6>c < 00 



r/2-4>o 



d^Si^-r-Oc)^ (53) 

COS(0) 



= -\ cos(arccos(^i /Oc) + M 



for 6*1 < 6*^ < 6*0 



for 6>i < 6>c < 6>o 



,(54) 



(55) 



sm(</>) 



for (p > 'k/2 — (po and 6 < C2 



(56) 



And: 



Plone: — 



Oi 



- 1 , for 6>2 < 6>c < 6>o (57) 



When the string "cuts the corner", entering and exiting via 
adjacent sides, Oc is: 



(58) 



Oo cos(0 - 200) - 20 
sin(20) ' 
for < 7r/2 — 00, ci < ^ < Co, 
and > 7r/2 — 0o, C2 < ^ < co 



Again using Eq. 48 we can find the component of Prad that 
comes from the string crossing the field via adjacent sides, P2: 

P,orner{Oc) = — d(j) (59) 

^ J^r^^r. 

) cos(0 — 0o) — 20 



^ /"^max 

2^ A... 



sin(20) 
sin(20) 



(60) 



0min and 0max mark the points when a second corner rotates across 
the string so that Oc reverts to either Oi or ^2. We can derive them 
from Eq. 50 and Eq. 51 and obtain: 



corner 



{Oc 



1 



sin(20), for6'c < Oi 



= — d(j) sin(20), 

^arccos(0i/0c) 

for Oi<Oc< O2 
= — dcj) sin(20), 

-'arccos(6»i /6»c) 

for 02<0c< Oo 

r(Oc) = -^,for(9c<^i 
27r^ 



(61) 

(62) 

(63) 
(64) 

(65) 



1 

2^ 



^2 



Finally, we add Pshort, -Piong and Pcomer to obtain: 



Prad(^c) — 



for6>c < Oi 



- 1 , for 6I2 < 00 < 6I0 (66) 



(67) 



, for 01 < 61c < 



1 



27^2 \\9cJ ^/W^t 

for e2<ec< \Jei + 0, 
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